function y=fun(k3,k,m,l1,l2,l3) 
m1=10.1*10^9+29.7*10^9*exp(-0.8)+122*10^9*exp(-1.6)+l1*exp(-1.6)+l2*exp(-1.6);
m21=3.29*10^9+10.1*10^9*exp(-(0.8+2/3*k3))+...
    3.29*10^9*exp(-(0.8+2/3*k))+29.7*10^9*exp(-(1.6+2/3*k3))+...
    l3*exp(-(0.8+2/3*k))+122*10^9*exp(-(2.4+2/3*k3))+...
    l3*exp(-2*(0.8+2/3*k))+29.7*10^9*exp(-(2.4+2/3*k3+2/3*k));
m22=l1*exp(-(2.4+2/3*k3))+122*10^9*exp(-(3.2+2/3*k3+2/3*k))+...
    29.7*10^9*exp(-(3.2+2/3*k3+4/3*k))+l3*exp(-3*(0.8+2/3*k));
y=17.86*k3*(1-exp(-(0.8+k3)*2/3))*m1/(0.8+k3)+...
    22.99*k*(1-exp(-(0.8+k)*2/3))*(m21+m22)/(0.8+k);
